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Abstract 

The persistent current I in integrable models of multichannel rings with both 
short- and long-ranged interactions is investigated. I is found to oscillate in 
sign and increase in magnitude with increasing interaction strength due to 
interaction-induced correlations in the currents contributed by different chan- 
nels. For sufficiently strong interactions, the contributions of all channels are 
found to add constructively, leading to a giant enhancement of I. Numer- 
ical results confirm that this parity-locking effect is robust with respect to 
subband mixing due to disorder. 
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The absence of macroscopic persistent currents in normal metal rings even at T = 
is a consequence of a parity effect stemming from Fermi statistics, first discussed by Byers 
and Yang fl[]: For N spinless electrons in a purely one-dimensional (ID) ring, the persistent 
current / is diamagnetic if N is even and paramagnetic if N is odd, independent of disorder 
and interactions @,@], and takes a maximum value Jo = evp/L for a clean ring, where Vp is 
the Fermi velocity and L the circumference of the ring. In a ring with many independent 
channels, I is the sum of many such diamagnetic and paramagnetic contributions, and is 
thus very small |T[]. Biittiker, Imry, and Landauer M argued that persistent currents should, 
nonetheless, exist at the mesoscopic level, and calculations with noninteracting electrons || 
predicted typical currents of order IqI/L for a diffusive metallic ring, where i -C L is the elas- 
tic mean free path. Subsequently, persistent currents were observed experimentally in both 
normal metal ||f|]] and semiconductor [|]J§ rings. Surprisingly, the persistent currents ob- 
served in metallic rings were roughly two orders of magnitude larger than those predicted by 
theories neglecting electron-electron interactions being of order 1$ in individual rings 

jTj. Mean- field calculations |TT|-[T3|1 for multichannel rings, renormalization group results 



for ID rings |L4|], and exact diagonalization studies |T|] found that the persistent current 
in disordered systems can be enhanced by repulsive interactions due to the suppression of 
charge fluctuations, but were unable to account for the experimentally observed magnitude 
of /. 

In this Letter, we propose a novel interaction effect which can lead to a large enhancement 
of / in multichannel rings, even in the ballistic regime, due to the fact that the parities 
of different channels are no longer independent in an interacting system. We investigate 
two integrable models of interacting M-channel rings, SU(M) fermions with inverse-square 
V(x) = g/x 2 and delta-function V(x) = U5(x) interactions. In the prior model, we show that 
for g > the persistent currents of all channels add constructively provided kpL > 2irM, 
leading to a large persistent current whose parity depends only on the total number of 
electrons N. In the SU(M) delta-function gas, / is found to oscillate in sign and increase in 
magnitude with increasing U due to a progressive condensation of electrons into the lowest 
subband, leading to a parity-locking effect for U > 16M A hvp/37T, where Vp is the Fermi 
velocity. A disordered two-channel ring with interchannel interactions is also investigated 
numerically, and shows, importantly, that the parity-locking effect persists even when the 
subbands are mixed strongly by disorder. Qualitatively, the parity-locking effect arises in a 
thin ring because for sufficiently strong repulsive interactions, electrons can no longer pass 
each other. The elimination of transverse nodes in the many-body wavefunction (which are 
necessary for two electrons to pass each other) leads to a state whose parity depends only 
on the total number of electrons, as discussed by Leggett 0, and whose persistent current 
is consequently large. 

Interacting spinless electrons in a non-disordered ring with M transverse channels, 
threaded by an Aharanov-Bohm flux (Hc/e)(j), may be represented by ID SU(M) fermions. 
The transverse degrees of freedom may be represented by an SU(M) spin variable a = 
1, . . . , M. In the absence of disorder, and for interactions which depend only on the elec- 
trons' coordinates along the ring (i.e., for thin rings), the number of electrons K a in each 
channel is conserved. The Hamiltonian of the system is 

i N a2 M 

^ = -oEtt2+E V{ Xi - X,) + £ K^, (1) 
z i=l ux i i<j cr=l 
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where N = J2a K a is the total number of electrons and e a is the energy minimum of subband 
a. Units with h — m — 1 are used. The Aharanov-Bohm flux leads to the twisted boundary 
condition [p]] 

^(xi<7i, ■■■ ,(xi + L)a h ■■■ , x N a N ) = e^{x 1 a 1 , • • • , XiO iy • • • , x N a N ). (2) 

For simplicity, let us consider equally spaced subbands e a+ \ — e a = A = EpjM. The 
subband splitting A plays the role of an SU (M) magnetic field. As we shall see, the effect 
of repulsive interactions is to renormalize this effective field, causing a condensation of 
electrons into the lowest subband. At T = 0, the equilibrium persistent current is given by 
7(0) = —(e/h)dE /d(p, where E ((j)) is the ground state energy of Eq. (]]]), subject to the 
boundary condition For simplicity, we consider the case <fi — ir/2 (1/4 flux quantum) in 
the following. 

Let us first consider a model with long-range interactions: V(x) = g/d(x) 2 , where d(x) = 
(L/ir)\ sin(7rx/L)| is the chord length along the ring. This model was introduced and solved 
by Sutherland |16| for the case M = 1 and = 0, and is integrable for g > 0. For g = 0, 
this model reduces to noninteracting electrons, and one finds 

W'l^K- 1 ^.. (3) 

(T=l 

For M 3> 1, this leads to the well-known || result |/| ~ M 1 ! 2 !^ due to the random parities 
of the different channels. For g > 0, on the other hand, the ground state is highly degenerate 
in the limit L — > oo in the absence of SU(M) symmetry breaking (A = 0) due to the strong 
repulsion of the potential at the origin, which prohibits particle exchange. In a finite ring, 
these states differ in energy by at most nhvp/L due to boundary effects; all electrons will 
thus be condensed into the lowest subband for A > irhvp/L, i.e., for kpL > 2ttM, which 
is satisfied provided the ring is sufficiently thin. The ground state of the system in this 
"ferromagnetic" state has the Jastrow product form 



= exp 




i Xn 

sin | -7i 



%AJ 1 t J 



sin ( — ^~ 7r ) (4) 



for < (f) < 7r, where a = if iV is odd and a = ir if iV is even. Here A = yg + 1/4 — 1/2. 
One readily verifies that \1/ is an eigenstate of Eq. (HD, has the correct symmetry, and obeys 
the twisted boundary condition This eigenstate is a positive vector when the particles 
are ordered, and is therefore the ground state of the system. The ground state energy is 
found to be 

p (i\ vr 2 (A+l) 2 iV(iV 2 -l) N (<P-aV 
The corresponding persistent current is 

<=7)7r/V 

W4) = ^ - i-^fMIo. (6) 
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The condensation of all electrons into the lowest subband caused by the strong repulsive 
interactions thus leads to an enhancement of the typical persistent current by a factor of 
M 1 / 2 in the ballistic regime, due to the suppression of the parity effect of Byers and Yang. 

Let us next consider a model with short-ranged interparticle interactions V(x) = U5(x). 
This system is integrable for arbitrary U, and the eigenenergies with twisted boundary 
conditions may be determined from a straightforward generalization of the nested Bethe 



ansatz of Sutherland |L7j to the case 7^ 0. The energy of the system may be expressed as 



N M 

E = Y. k V 2 + H K ° £ °i ( 7 ) 

where the pseudomomenta kj are a set of N distinct numbers which satisfy the coupled 
equations 



N-Ki u. _ ^(1) +HJ/2 

exp[i(L*,-0)] = n I % r/ t ( § ) 
= i kj — Aa — iU/2 



N , A (n) _ A (n-1) _ lU/2 AWi A (n) _ A (n+1) _ Nn A M _ A^) - lV 

M A<?> - A?" 1} + iU/2 M - A^ +1) + iU/2 " ~ }k ^ - If + iU ' n - 2 '"--' M 

(9) 

where A^, a = 1, . . . , N n = N — J2a=i are distinct numbers, with A^ = kj. For A = 0, 
the ground state is an SU(M) singlet when is an odd multiple of M. As A increases, 
electrons in the higher subbands are transferred to lower subbands, until all electrons are 
condensed into the lowest subband for A > A c . This phenomenon is analogous to the spin- 
polarization transition of the ID Hubbard model in a magnetic field, studied by Carmelo et 
al. and by Frahm and Korepin [15 ] . Using the techniques of Ref . [IS ] , one finds 



A c = (1/4tt)[[/ 2 + (2vm) 2 ] tan^^vm/fT) - Un/2, (10) 

where n = N/L. For A > A c , the system is in a parity-locked state, with persistent current 
given by Eq. (|6|). It is useful to consider some limiting cases of Eq. flip]) : Using kp — im/M, 
one finds A c ~ M 2 E F for U = and A c ~ 1QM 3 E f v f /3ttU for > M. For fixed A, 

the critical interaction strength required to enforce parity-locking is thus 

U c ~ 16M 4 w f /3tt. (11) 

For M 3> 1, very strong interactions are thus required to cause complete parity-locking, 
which would lead to a macroscopic persistent current. This result is in contrast to that for 
the preceding model, which exhibited parity- locking for any value of the interaction g > 0, 
provided the ring was sufficiently thin. 

In order to see how the parity-locking effect develops as a function of interaction strength, 
let us consider the simple case of a two-channel ring with odd; then the parities of the 
two channels are necessarilly opposite. For mesoscopic rings, SU{2) excitations of the type 



considered by Kusmartsev [0 and by Yu and Fowler [ 20| , which lead to a 4>o/N periodicity 
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of the persistent current, can be neglected [20]. The persistent current for M = 2 is then 
given by 



I(0o/4) = (-l) Xl (^i " K 2 )I /N. (12) 

For A <C A c , the polarization (Ki — K 2 )/L ~ 2%A, where the susceptibility x may be 
evaluated from Eqs. using the method of Shiba fl21| ; one obtains x — 2/ttvp for U = 



and x — 3U/(2ttvf) 2 for U ^> uj?. The magnitude of the persistent current is thus given by 
|/| — ieA/nhN for U = 0, and is increased by the factor 3U/8nvp for Vp <C U < [4. As 
[/ is increased, the persistent current thus oscillates in sign and grows roughly linearly in 
magnitude due to the progressive transfer of electrons from the upper to the lower subband. 
For M 3> 1, the evolution of the system toward the parity- locked state as U is increased 
will of course be more complicated, but one nonetheless expects I to fluctuate in sign and 
increase in magnitude as electrons condense into the lowest subband. 

A peculiarity of the integrable models considered to this point is that the number of 
electrons in each channel is a constant of the motion. Both disorder and more realistic 
interactions which depend on the transverse coordinate will break this symmetry, and it is 
therefore important to verify that the parity-locking effect is not destroyed. To this end, we 
have considered a disordered two-channel ring, modeled in the tight-binding approximation, 
with a nearest-neighbor interchain interaction V included to induce interchannel correlations. 
The Hamiltonian is 



i=i 



2 ^ 

^ /L c\ a c i+la + H.c. + s ia p ia ) + — (4^2 + H.c.) + Vpiip i2 



(13) 



where c\ a creates a spinless electron at site i of chain a, p ia = c\ a c ia , and e ia is a random 
number in the interval [—e/2, e/2]. The interchain hopping determines the subband splitting 
A. Fig. 1 shows the persistent current for rings with 5 spinless electrons on 18 sites as a 
function of V, calculated using the Lanczos technique. Both the ensemble average of / 
for 500 systems (squares) and the persistent current of a typical system (solid curve) are 
indicated. The error bars indicate the standard deviation SI = ((I 2 ) — (J) 2 ) 1 / 2 over the 
ensemble (not the statistical uncertainty in the mean value). The subband splitting A = 0.8 
is chosen so that in the absence of disorder and interactions, K\ = 3 and K 2 = 2, leading 
to a large cancellation of the persistent current due to the different parities of the two 
channels. The on-site disorder e = 2 > A, Ep mixes the two channels, but does not lead 
to strong backscattering. For V = 0, the sign of I is essentially random, and (/) is slightly 
diamagnetic. As V is increased, (J) oscillates in sign and increases in magnitude, becoming 
strongly diamagnetic for large V. |(/)| is increased by a factor of 10 as V is increased from 



to 20, while y (I 2 ) is increased by a factor of 3. How can we understand this behavior, so 
analogous to that of the ballistic system with delta-function interactions described above, 
given the strong intersubband scattering? While the subband occupancies are no longer 
constants of the motion, there is a corresponding topological invariant in the disordered 
system, namely, the number of transverse nodes in the many-body wavefunction (i.e., 
nodes which encircle the AB flux <fr). The lowest subband has no such nodes, while each 
electron in the second subband contributes one transverse node. In order for two electrons 
to pass each other as they circle the ring, such a transverse node must be present. As 
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V increases, it becomes energetically unfavorable for electrons to approach each other, so 
transverse nodes in the many-body wavefunction will tend to be suppressed. In the strongly- 
correlated limit, all such nodes will be eliminated, leading to a state whose parity depends 
only on the total number of electrons. In such a state, the persistent currents of all channels 
add constructively, leading to a large persistent current (see Fig. 1). It should be emphasized 
that while the enhancement of the persistent current shown in Fig. 1 is relatively modest, a 
much larger enhancement would be expected in a system with M 3> 1, based on the above 
arguments. 

In conclusion, we have proposed a novel interaction effect which leads to a large enhance- 
ment of the persistent current in multichannel rings due to correlations in the contributions 
of different channels. It was shown that even when interactions are weak compared to those 
necessary to enforce complete parity-locking, as is likely to be the case in metallic rings such 
as those studied in Refs. the persistent current may still be substantially enhanced 

by interchannel correlations. It should be emphasized that the parity-locking effect holds 
in both ballistic and disordered systems; it is therefore complementary to mechanisms pre- 
viously proposed for the enhancement of the persistent current [pTT]— 115[1 , which rely on the 
competition between disorder and interactions. It is likely that both such mechanisms are 
important to explain the anomalously large observed value |]J7|] of the persistent current in 
disordered metallic rings. 

C. A. S. thanks M. Biittiker and F. Hekking for valuable discussions. This work was 
supported by the Swiss National Science Foundation. 
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FIGURES 



FIG. 1. Persistent current / = — (e/ft)<9£7o/d0|</>=0 o /4 of a disordered two-channel ring with 5 
spinless electrons on 18 sites as a function of the interchain interaction V. The current is given 
in units of Iq = evp/L. The amplitude of the on-site disorder is e = 2 and the subband splitting 
is A = 0.8. Solid curve: persistent current for one realization of disorder. Squares: ensemble 
average (I) for 500 systems. The error bars indicate the width 51 = ({I 2 ) — (I) 2 ) 1 / 2 of the current 
distribution. Note that the persistent current is diamagnetic for large V, as expected for a system 
with N odd due to parity-locking; this is true for all realizations of disorder. 
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